Vortex penetration and flux relaxation with arbitrary initial conditions in non-ideal 

and ideal superconductors 
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Vortex penetration and flux relaxation phenomenon carry the information about the pinning abil- 
ity, and consequently current-carrying ability, of a type-II superconductor. However, the theoretical 
descriptions to these phenomena are currently limited to the cases with special initial conditions. 
A generalization to the recently developed infinite series models is presented here. It is shown that 
one can convert a vortex penetration process with a non-zero initial internal field into a process with 
a zero initial internal field by introducing some time parameters. Similarly, one can also convert a 
flux relaxation process starting with an arbitrary internal field into a process starting with a melting 
internal field by introducing a virtual time interval. Therefore, one can predict the melting internal 
field (or critical current density) from a flux relaxation process starting with a lower internal field. 
Finally, it is shown that the vortex penetration process in an ideal superconductor is strongly time 
dependent because of the surface barrier and internal field repulsive force. But the flux relaxation 
process does not occur in the ideal superconductor. 

PACS numbers: 74.25.Op, 74.25.Uv, 74.25.Wx 



I. INTRODUCTION 

In vortex penetration process, the vortices are pushed 
into a type-II superconductor by external driving 
forced—. The internal field is then an increasing func- 
tion of time. At lower applied fields or lower tempera- 
tures, the vortex penetration can be regarded as a process 
of vortices hopping between adjacent pinning centers** 
which can be described by the Arrhenius relation. By 
proposing an internal field dependent activation energy, 
one can obtain the corresponding time evolution equation 
of the internal field*** In application, a superconducting 
device may be repeatedly loaded and unloaded, and part 
of the internal field remains in the bulk of the super- 
conductor because of pinning. 12 In this case, the vortex 
penetration process starts with a nonzero initial inter- 
nal field. However, the existing theoretical models** are 
only applicable to the vortex penetration process with 
zero initial internal field. Thus, we need to extend these 
theoretical models and make them applicable to the vor- 
tex penetration processes with a nonzero initial internal 
field. 

On the other hand, in the flux relaxation process, the 
internal field (or current density) of a superconductor 
decays due to the spontaneously vortex hopping between 
adjacent pinning centers »**~— The internal field (or cur- 
rent density) is then a decreasing function of time. Using 
the time evolution equation of the internal field, one can 
determine the melting internal field from a flux relax- 
ation measurement by applying an external field above 
the melting value. But this measurement may be tech- 
nically difficult at lower temperatures where the melting 
internal field is large. 1 - Therefore, we need to find a gen- 
eral formula which can predict the melting internal field 
from a flux relaxation process starting with a lower initial 
internal field. 

Furthermore, it is known that ideal superconductors 



(clean superconductors)*^ have a perfect crystal lattice 
and are free of pinning centers. This indicates that mag- 
netic vortices cannot be pinned down in an ideal super- 
conductor. But it is not clear how fast the vortices can 
penetrate into or hop out of the ideal superconductor. Is 
the time evolution equation of an internal field logarith- 
mic or non-logarithmic? Thus, it is desirable to develop 
a mathematical model to describe the time dependence 
of an internal field in the vortex penetration and flux 
relaxation process of an ideal superconductor. 



In this article we considered the fact that, in the flux 
relaxation and vortex penetration process, the pinning 
potential of the vortex does not change but the activa- 
tion energy doesi*^ Therefore, we constructed the general 
equations for time dependence of activation energy and 
of the internal field by introducing some equivalence time 
parameters related to the pinning potential. Using these 
equations, one can predict the maximum internal field 
in a vortex penetration process and predict the melting 
internal field from a flux relaxation process with an arbi- 
trary initial internal field. We also showed that the vortex 
penetration into an ideal superconductor is strongly time 
dependent. 



Because the activation energy is a function of coher- 
ence length and penetration depth, it includes informa- 
tion about the anisotropy of a material. Therefore, we 
intend to construct the theoretical models which can be 
applied to both high-T c and low-T c superconductors. In 
case that the anisotropy of the material has to be spec- 
ified, we discuss it separately (Section IV). Also, in the 
derivation we only considered bulk pinning potential and 
Bean-Livingston surface barrier, but ignored geometry 
barriers. 
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II. VORTEX PENETRATION PROCESS WITH 
ARBITRARY INITIAL CONDITIONS 

The purpose of this section is to find a general expres- 
sion for the time dependence of an internal field in a vor- 
tex penetration process with an arbitrary initial internal 
field. To do this, we need the internal field dependence 
of activation energy and time dependence of activation 
energy. 
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A. Internal field dependence of activation energy 
in vortex penetration 

The Bean model^i shows that, under an applied field 
B a , the internal field should approach a maximum value 
B e when the superconductor reaches an equilibrium 
state. Because of the surface screening effect (or Meissner 
effect), B e must be smaller than B a , that is, B e < B a . 

a. Series expression of activation energy In vortex 
penetration process, the repulsive force of the internal 
field B prevents vortex motion and reduces vortex hop- 
ping rate. The activation energy U p is then an increasing 
function of B. Thus, dU p /dB > 0. In early studyii, U p is 
expressed as a series of B (or magnetization M). For the 
convenience of calculating the equilibrium field B e , here 
I expand U p as a series of the normalized field a = B / ' B e , 
that is, 
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U p (a) = U p0 + y2aia l , 



(1) 



Time 

FIG. 1: Schematic diagram of U p , the activation energy of 
a vortex penetration process. U c is pinning potential, U p o is 
the activation energy at vanishing internal field, Ui is initial 
activation energy, U e is equilibrium activation energy, Bi is 
initial internal field and B e is equilibrium internal field, (a) 
Internal field dependence of activation energy. U p is an in- 
creasing function of internal field B. (b) Time dependence of 
activation energy. U p is an increasing function of time t. 



where U p q = Ubl + U c is the activation energy at van- 
ishing internal field, Ubl is the Bean-Livingston surface 
barrier,— and U c is pinning potential. 

b. Activation energy at initial internal field If the 
initial internal field Bi > (at t = 0), then the cor- 
responding initial activation energy Ui > U p q. Using 
Eq.([T]), we can express Ui in terms of <7j = Bi/B e , that 
is, 



Ui = u 



pO 



i=i 



(2) 



c. Activation energy at equilibrium internal field At 
the equilibrium field B e , the corresponding equilibrium 
activation energy is 



U e = U p0 + 



1=1 



(3) 



On the basis of the above discussion, we can now draw 
a schematic diagram of U P (B), the internal field depen- 
dence of activation energy in a vortex penetration pro- 
cess, as shown in Fig. 1(a). 

Inverting Eq.(fT|) and using definition a = B/B ei we 
obtain the expression of internal field B in terms of the 



activation energy U p , that is, 

oo 

B{U p ) = B e J2h{U p -U p0 ) 1 , 
1=1 

where the coefficients 6/ areii 
1 1 



(4) 



l(l + l)---(l-l + s + t + u + ---) 



sltlul ■ ■ ■ 
aij \aij \ai 



(5) 



and s + 2t+3u + - ■ ■ = I — I. On considering the symmetry 
between Eq.((T]) and Eq.Q, we can obtain the inverse 
coefficients a/ by doing a commutation to the coefficients 
bi o a/. 

Eq.([4]) shows that, to obtain the time dependence of 
internal field B(t), we still need the time dependence of 
activation energy U p (t). Let us discuss this in the next 
section. 
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B. Time dependence of activation energy in vortex 
penetration 

1. Solution of the time dependence of activation energy 

We have shown that the activation energy of the vortex 
penetration process, U p , is an increasing function of the 
internal field B. Because B is an increasing function time 
t (dB/dt > 0), Up is an increasing function of t. Using 
the Arrhenius relation, we have^i 



dB 
~dt 



Ce -U P /kT 



(6) 



where C is a positive proportional constant, k is the 
Boltzmann constant and T is temperature^ 

Consider a superconductor with an initial internal field 
Bi and initial activation energy fjj. Now apply an exter- 
nal field B a to the superconductor at zero time (t = 0). 
As the time increases to t, the activation energy increases 
from Ui to Up, Rewrite Eq.© and integrate it on both 
sides 



Ui 



dB 



With logarithmic accuracy, we obtain the following 
equation^ 



U p (t) = kTln e 



T 



(7) 



where r = kT/[C(dU p /dB)] is a short time scale 
parameter—. Eq.(J7|) describes the general time depen- 
dence of activation energy in a vortex penetration pro- 
cess, U p (t), which is an increasing function of time t. 



2. Meaning of virtual time intervals 

To assist physical understanding and easy calculation, 
let us define the following time parameters 



ti = T 



i/kT 



1 



tv = ti - t a = t (e 



Ui/kT _ M p0 /kT 



(8a) 
(8b) 
(8c) 



where U p o is the activation energy at vanishing inter- 
nal field (Eq.flTJ)) and f/, is the initial activation energy 
(Eq.©). Thus, Eq.flTJ) becomes 

U p (t) = kTln (l + *±*) = kTln (l + to + t J +t 

(9) 

a. Meaning of t v In Eq.©, if we let t — > 0, then 
U p Ui (use Eq.flHbJ) and Eq.flgcJ)). If we further let 
t v — > 0, then Ui — > U p o- This indicates that the physical 



meaning of t v is: the virtual time interval during which 
the activation energy U p (t) increase from U p q to Ui. 
Putting t v = in Eq.©, we have 



U p (t) = kTln 1 + 



ta + t 



(10) 



Eq. ([8c|) shows that, if t v — 0, then Ui — U p0 - From 
Eq.([2]), we have the initial internal field Bi = 0. Thus, 
Eq. lfTU)) is the time dependence of activation energy in 
a vortex penetration process with zero initial internal 
field.il 

By doing a time transformation t' = t v + t, we can 
convert Eq.Q into the form of Eq. (fT0| . In other words, 
we can convert a vortex penetration process with nonzero 
initial internal field into a process with zero initial inter- 
nal field by introducing a time parameter t v as defined in 
Eq. (l5c| . This can be clearly seen in the later discussion 
on the time dependence of internal field. 

b. Meaning of to Because U p o is the activation en- 
ergy at vanishing internal field (Eq.([T])), a superconduc- 
tor with U p o = means that it is free of pinning. The 
superconductor is then an ideal superconductor, or clean 
superconductors^ (Here I assume that the superconduc- 
tor is isotropic. The layered superconductor is discussed 
in section IV). Thus, using Eq. (l8bl) we can conclude that 
the physical meaning of to is: the time parameter equiv- 
alent to the potential U p o in a non-ideal superconductor. 

Putting t = in Eq. pU)) . we have 



U p (t) = kTln 1 



t 



(11) 



t v = means that the initial internal field Bi = 
(Eq. flTUl) '). and t — means that the superconductor 
is an ideal superconductor (Eq.([8b|). Thus, Eq. flTTj) is 
the time dependence of the activation energy of a vortex 
penetration process in an ideal superconductor with zero 
initial internal field. 

By doing a time transformation t" = to + t, we can 
convert Eq. (fT0l) into the form of Eq. (fTTj) . In other words, 
we can convert a vortex penetration process of a nonidcal 
superconductor into a process of an ideal superconductor 
by introducing a time parameter to, as defined in Eq. (|8bl) . 

From the above discussions we know that there is a 
equivalence relation between the time t and activation 
energy U p . To see this, let us invert Eq.flTTj), 



t = T 



f e U P /kT _ ^ 



(12) 



Eq. (|12p shows that, in an ideal superconductor with 
zero initial internal field, the activation energy of a vortex 
penetration process can be mapped into a time parame- 
ter, that is, f :U p ^t,Up£ [0, U e ]. Although Eq.flJJJ) is 
obtained from an ideal superconductor, it can generate 
all the relations for nonideal superconductors because the 
potential U p q is equivalent to a time parameter to. 

The time dependence of activation energy in a vortex 
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penetration process is schematically shown in Fig. 1(b). 



C. Time dependence of internal field in vortex 
penetration 

The time dependence of internal field B{t) can be ob- 
tained by combing the internal field dependence of activa- 
tion energy U p (a) with the time dependence of activation 
energy U p (t), 



1. Infinite series activation energy 
Substituting Eq.® into Eq.flJ), we have 

oo 

B(t) = B^Hw^t)} 1 , 

where 



1=1 



w p (t) = kTln ( 1 + ° — v ^ I — Upo- 



(13) 



(14) 



Eq.(fT3|) describes the time dependence of the internal 
field B(t). Fitting Eq. flU]) to the experimental data, we 
can obtain various parameters and then calculate the ac- 
tivation energy^ 

It should be mentioned that, to determine B e , we need 
both Ea.([i"3]l and the inverse of Eq.©, 



which the activation energy increase from C/ p o to the ini- 
tial value Ui. Now we can explain it in terms of internal 
field: t v is the time interval during which the internal 
field increases from to the initial value Bi. 

Eq. (fT3l also shows that a vortex penetration process 
with an initial internal field Bi can be converted into a 
process with zero initial internal field by introducing a 
virtual time interval t v . 

In Eq. (j8"a|) and Eq. ([5c|) we express the time parameters 
ti and t v in terms of the initial activation energy Ui. But 
in practice the measurable physical quantity is the initial 
internal field Bi. Therefore, it should be more convenient 
to express ti and t v in terms of Bi (o~i — Bi/B e ). To do 
this, a direct way is to invert Eq. (|17[) . but an easier way 
is to use Eq.fl2J|. 

Substituting Eq.® into Eq. (|8a| . we have 



,U p0 /kT 



i=i 



(18) 



Substituting Eq.© into Eq.([8c|). we have 



_ Tp U p0 /kT 



t v = re 



11 = 1 



(19) 



If we further put to = in Eq. ((T7|) . we still have Bi = 
0. Therefore, in a nonideal superconductor, the physical 
meaning of to cannot be explained in terms of an internal 
field. 



l = ^bi(U e -Uj 



p0) 



(15) 



i=i 



This can be obtained by simply putting B — B e and 
U p = U e in Eq.®. 

If we are only interested in finding out the time de- 
pendence of the internal field, then we can simplify 
Eq. fllUl) by absorbing B e into the coefficients b\. Putting 
hi = biB e in Eq.flT^f. we have 



B(t)=J2h l [w p (t)} 1 . 



(16) 



i=i 



This is equal to expanding the activation energy as a 
function of the internal field B (Ref. [Til) . that is, U P (B) = 
Upv+YZiaiB 1 . 

Putting t — in Eq. (fT"3| . we obtain the initial internal 
field 



Bi 



B e J2 b i 



1=1 



kTln 1 



to + tv 



U, 



/jo 



(17) 



Putting t v = in Eq. p7| and using Eq. (|8bp . we have 
Bi =0. This is consistent with the discussion on Eq. flTUl) . 
in which we interpreted the physical meaning of t v using 
the activation energy, that is, the time interval during 



2. Quadratic activation energy 

Without inelastic deformation^, a noninteracting elas- 
tic vortex system has ai = (I > 2). Thus, Eq.(fT|) re- 
duces to the following quadratic activation energy 



U P (B) = Upo + a x o + a 2 o- 2 , 



(20) 



and Eq.Q reduces to U e = U p o + ai + a 2 . 

Substituting Eq.© into Eq. (|2"0|) . we have (choose one 
of the solutions that is an increasing function of time) 



W 2a 2 



1 + 4— w p {t) - 1 



(21) 



where w p (t) is defined by Eq. (fl"4"l) . 
Putting t = in Eq.([2T|l. we have 



Bi — B, 



ai 
'2a 2 




UTln ! 1 + 1 - U p0 



Eg. ([18]) reduces to 

£. — T e (U p0 +ai(Ji+a2<j'l )/kT _ ^ 



- 1 

(22) 
(23) 
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and Eq. (fT9|) reduces to 



t v = re 



U p0 /kT 



^(a 1 a i +a 2 a'f)/kT 



(24) 



3. Linear activation energy 



In a noninteracting rigid vortex system, the elastic de- 
formation vanishes^. Thus, the coefficient ci2 = and 
Eq.([T]) reduces to the following linear activation energy 



U P (B) = U pQ + aid, 

and Eq.((3|) reduces to U e = U p o + a\. 

Substituting Eq.© into Eq.(|2l)]). we have 

B(t) = — u)p(t). 

where w p (t) is defined by Eq. ([In- 
putting t = in Eq.(|26|). we have 



(25) 



(26) 



B, 



Be 



kTln 1 



to + t, 



U, 



pO 



Ea.(fl8|) reduces to 



U = t 



and Eq. p^|) reduces to 

t v = T e u -°' kT (e a ^' kT - I 



(27) 



(28) 



(29) 



III. FLUX RELAXATION PROCESS WITH 
ARBITRARY INITIAL CONDITIONS 

In the study of flux relaxation, a vortex's activation en- 
ergy is usually expressed as a function of current density 
and the corresponding time dependence of current den- 
sity can be derived. 17 ' 23 ~ — In the present work, we are 
going to discuss the flux relaxation process in an ideal 
superconductor. Because the flux relaxation process is 
essentially a process of vortex motion, it should make 
more sense to use internal field instead of current den- 
sity. In case of need, one can convert the formulas in 
terms of internal field into formulas in terms of current 
density or magnetization by doing a replacement. 



A. Internal field dependence of activation energy 
in flux relaxation 

At a temperature T, a vortex system will melt as the 
internal field increase to a value B m . Corresponding to 
this B m , there should be a melting field as defined in the 
conventional magnetic phase of type-II superconductors. 



The melting field is not an internal field, but an external 
field. 

It should be emphasized that the melting internal field 
B m is different from the equilibrium field B e in the vor- 
tex penetration process as discussed in Eq. ([I} . In the flux 
relaxation process, B m is the maximum internal field for 
the vortex system to keep a lattice structure^ (corre- 
sponding to the external melting field), which is deter- 
mined by the temperature, pinning and material of the 
superconductor. But in a vortex penetration process, B e 
is the maximum magnetic field can penetrate into the su- 
perconductor under an applied magnetic field B a , which 
is determined by the applied field B a . 

a. Series expression of activation energy In the flux 
relaxation process, the repulsive force between vortices 
can enhance the vortex motion and increases vortex hop- 
ping rate. The activation energy U r is then a decreas- 
ing function of the internal field B. This indicates that 
dU r /dB < 0. According to an early study-i 7 -, we can 
express U r as a series of internal field B. For the con- 
venience of calculating the melting internal field B m , I 
expand U r as a series of the normalized internal field 
A = B/B m , that is, 



U r (X) = U r0 



OO 

E 

1=1 



ci\ l 



(30) 



where U r o — U c — Ui m . The parameter U c is the pinning 
potential inside the bulk and Ui m = (<&q/A-kX) 2 K (2x/X) 
is a reduction to the activation energy caused by the sur- 
face imaging force. 17 ' 22 Eq. (|37))) shows that U r o is the 
activation energy of vortices at vanishing driving force 
(B = 0) and is also the maximum activation energy of 
the entire flux relaxation process. 

b. Activation energy at initial internal field Let Bi 
be the initial internal field at time t = and Ui be the 
corresponding initial activation energy. Putting B = B- t 
(or A<j = Bi/B m ) in Eg. (|3"0"1) , we have 



Ui = U, 



r0 



E 



(31) 



Since the activation energy U r is a decreasing function 
of the internal field, Eq. (|3"Tj) indicates that U < U r o- 

c. Activation energy at melting internal field At the 
melting internal field B m , the activation energy U r is 
zero. Eg. (f3"0]) gives 



U r Q 



1=1 



(32) 



This is a constraint condition for the coefficients c/. 

On the basis of the above discussions, we can now draw 
a schematic diagram of U r (B), the field dependence of 
activation energy in a flux relaxation process, as shown 
in Fig. 2(a). 

Inverting Eg. ([3H)l and using the definition A = B/B m , 
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Time 



FIG. 2: Schematic diagram of U r , the activation energy of 
a flux relaxation process. U c is pinning potential, U r o is the 
activation energy at vanishing internal field, Ui is initial ac- 
tivation energy, Bi is initial internal field and B m is melting 
internal field, (a) Internal field dependence of activation en- 
ergy. U r is a decreasing function of internal field B. (b) Time 
dependence of activation energy. U r is an increasing function 
of time t. 



we obtain the expression of internal field B in terms of 
activation energy U r , that is, 

oo 

B(U r ) = B m di(U ra - U r ) 1 , (33) 
1=1 

where the coefficients di and cj has the same relation as 
that of b t and a/ (See Eq.©). 

Eq. (j3"3"|) shows that we can obtain the time dependence 
of the internal field B(t) by finding out the time depen- 
dence of the activation energy U r (t). This is further stud- 
ied in the next section. 



B. Time dependence of activation energy in flux 
relaxation 



where C is a positive proportional constant. 

Consider a superconductor in which an arbitrary inter- 
nal field starts to decay at zero time (t = 0). The initial 
activation energy of the vortices is Ui. As time increases 
to t, the activation energy increases to U r . Similar to 
Eq. © , we obtain the following equation^ 

U r (t) = kTln(l + ^^j , (35) 

where 

U = r{e^l kT -l) (36) 

and t — —kT/[C(dU r /dB)]. Here we include a negative 
sign to ensure that r is a positive number because in 
the flux relaxation process the activation energy U r is a 
decreasing function of the internal field B, i.e., dU r /dB < 
0. 

Eq.([35| shows that the activation energy of a flux re- 
laxation process starting with an internal field below the 
melting internal field can be converted into that of a pro- 
cess starting with the melting internal field by introduc- 
ing a virtual time interval U. The physical meaning of t% 
is: the virtual time interval during which the activation 
energy U r (t) increases from to the initial value Ui. 

It should be mentioned that the time parameter tj 
in Eq . (|35|) can not be split into to and t v as we did in 
Eq.©. This is because in the flux relaxation process, 
the activation energy at vanishing internal field, U r o, is 
the maximum activation energy. The initial activation 
energy Ui satisfies Ui < U r Q. Thus, Ui cannot be further 
split. But in vortex penetration process the activation 
energy at vanishing internal field, U p o, is the minimum 
activation energy. The initial activation energy Ui satis- 
fies Ui > U P Q. Thus, Ui can be split into two parts. 

Eq. (|3"5)) describes the general time dependence of the 
activation energy in a flux relaxation process, U r (t), 
which is an increasing function of time t. It turns out 
that both the activation energy of the flux relaxation pro- 
cess U r (t) (Eq. (f3"5|) ') and the activation energy of vortex 
penetration process U p (t) (Eq.©)are increasing function 
of time t. 

The time dependence of activation energy in a flux 
relaxation process is schematically shown in Fig. 2(b). 



In the flux relaxation process, the vortex's activation 
energy U r is a decreasing function of the internal field B. 
Because B is a decreasing function of time t (dB jdt < 0) , 
U r is an increasing function of t. Similar to Eq.©, we 
have 

^ = -Ce~ u ^ kT (34) 



C. Time dependence of internal field in flux 
relaxation 

The time dependence of the internal field B(t) can be 
obtained by combing the time dependence of activation 
energy U r (t) fEq. (|3"5")) ) with a detailed internal field de- 
pendence of activation energy U r (B). 
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1. Infinite series activation energy 



2. Quadratic activation energy 



The infinite series activation energy^ is shown in 
Eq. (f3"U| and its inverse is shown in Eq. ([55]) . Substituting 
Eq.([3"5"j) into Eq. (l3"3"l) . we obtain the time dependence of 
the internal field 



B{t) = B m J2Mwr(t)} 1 , 



i=i 



where 



w r {t) = U r o — kTln 1 + 



U + t 



(37) 



(38) 



Putting t = in Eq. (|3"Tl) . we obtain the initial internal 
field 



Bi — B m ^ di 



U r0 - kTln 1 + 



u 



(39) 



Putting ti = in Eq. (|36|) . we have Ui = 0. It means 
that the flux relaxation process starts with the melting 
internal field, i.e., Bi = B m . Thus, Eq. (|3"9")l gives 



(40) 



In a noninteracting elastic vortex system with vanish- 
ing inelastic deformation—^, the coefficient q =0 (I > 2). 
Thus, Eq.(l3U)) reduces to the following quadratic activa- 
tion energy 



U r {X) = U ra - ciA - c 2 A 2 , 
where c\ + c 2 = U r o (Using Eq. 



(42) 



Substituting Eq. (l3"5)l into Eq. (|4^|) . we obtain the time 
dependence of internal field B(t) (choose one of the so- 
lutions which is a decreasing function of t) 



B{t) = B v 



2c 2 



, °2 



1 + A^jw r {t) - 1 



where the function w r (t) is defined in Eq.([38 
Putting t = in Eq.g2]), we have 



Bi — B 7 



Cl 

2c 2 



and Eq. ([4T|) reduces to 



L = T 



1 = 1 



,(U r0 -ci\i-c 2 \?)/kT _ 1 



(43) 





U r0 - kTln ^1 + ^-J 


- 1 









(44) 
(45) 



Eq. pO)) is a constraint condition of the coefficients di. 
It is also the inverse of Eq. (l3"2l (rewriting as U r o — 
Ei=iQ x 1). 

Eq. ([3T|) describes the time dependence of the internal 
field B(t). It shows that a flux relaxation process starting 
with an initial internal field Bi can be converted into a 
process starting with the melting internal field B m by 
introducing a virtual time interval t<. Fitting Eq. (|3T|) to 
experimental data, we can obtain the melting internal 
field B rn . This means that we can calculate B m from a 
flux relaxation process with an arbitrary initial internal 
field Bi. 



In Eq. (l36l) . we expressed the virtual time interval ti 
in terms of the initial activation energy Ui. Let us now 
express ti in terms of the initial internal field Bi (Aj = 
Bi/B m ). Substituting Eq.(|3l)l into Eq.([36j). we have 



U = t 



e U r0 /kT TT e -ciX\/kT _ 1 
1=1 



(41) 



When deriving Eq. (|3"5|) . we explained the physical 
meaning of the virtual time interval ti using activation 
energy. Now, we can explain it using the internal field 
from Eq. (|37[) : U is the time interval during which the 
internal field reduces from the melting value B m to the 
initial value Bi. 



3. Linear activation energy 



In a noninteracting rigid vortex system^, the coeffi- 
cient c 2 = 0. Eq. (|4"2")l reduces to the following linear 



activation energy^ 

f7 r (A) = U rQ - ciA, 

where a = U r o (Using Eq. ([32|) ). 

Substituting Eq. (|33|) into Eq. pp]) . we have 

B(t) = —W r {t). 
Cl 

where the function w r (t) is defined in Eq. 
Putting t = in Eq.gT]), we have 



Bi = 



B m 

Cl 



r,„ - kTln ( 1 + - 



and Eq. (|4"Tj) reduces to 



(46) 



(47) 



(48) 



(49) 
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IV. DISCUSSION 



A. Vortex motion in ideal superconductors 



An ideal superconductor has perfect crystal lattice and 
is free of defects^ 9 But for a layered superconductor, even 
if it is "ideal" , the space between the layers can still be- 
have as pinning centers^ Thus, a vortex can be pinned 
down when it lies between two layers. In this section I 
am only interested in studying a physical system in which 
the vortices are free to move. Therefore, I assumed that 
the ideal superconductor under consideration is isotropic. 
It is not only free of pinning centers, but also free of any 
pinning potential, i.e., U c = 0. 



Vortex penetration 



In an ideal superconductor, the pinning potential is 
U c = 0. Thus, the internal field dependence of activation 
energy of a vortex penetration process, Eq.([T|), reduces 
to 



CO 

U p (a) = Ubl + ^cicf 1 . 

i=i 



(50) 



Eq. (|50[) is normal, which indicates that the vortex pen- 
etration process may occur in an ideal superconductor. 

Let us now check the time dependence of an internal 
field in the ideal superconductor. First, use Eq.([7]) and 
rewrite Eq. flTBI as 



,(Ui-U p0 )/kT 



i=i ^ 

Putting U c = in Eq.(f5"Tj). we have 



+ 



t 



?U p o/kT 



(51) 



B(t) = ^hi\ kTl 



e {Ui-U BL )/kT + 



t 



,U B L/kT 



(52) 

Since {U l -U BL ) > 0, we have e (^-^i,)/fcT > L Thnfj ^ 
Bit) > 0. The simulation of Eq. ([52|) is shown in Fig. 3. 

One can see that, in an ideal superconductor, the in- 
ternal field B is an increasing function of time t. This 
shows that although an ideal superconductor does not 
have pinning centers, the vortex penetration into this su- 
perconductors is strongly time dependent. The reason is 
that the repulsive force of the internal field (a field gradi- 
ent formed in the bulk) prevents the vortices penetrating 
into the bulk of the superconductor. The vortex motion 
is then retarded. 



o 



Ideal Superconductor 




Simulation of 
Vortex Penetration 



10~ 1 10° 10 1 10 2 10 3 
Time (s) 



10" 



FIG. 3: (Color online) Simulation of vortex penetration into 
an ideal superconductor. The time dependence of an internal 
field is simulated with equation: B p (t) = h\w p (t) + fetOp(t), 
where w p (t) = (25k) x ln(l + t), hi = 8 X 10~ 5 /(25A:), h 2 = 
-4.5 x 10" 6 /(25fc) 2 (See Ea.(T52)l). Here I choose a position 
away form the surface, where Ubl = 0. Also, I choose r = 1 
and Bi = 0. In an ideal superconductor the pinning potential 
is U c — 0, we have U = Ubl = 0. 



2. Flux relaxation 



In an ideal superconductor, the pinning potential is 
U c = 0. Thus, the field dependence of the activation 
energy of a flux relaxation process, Eq. (|3"0T) , reduces to 



Ur(X) 



-U im ~^2ci\ l 
1=1 



(53) 



This shows that U r (X) < 0. However, for a flux relaxation 
process to occur, the activation energy must be positive. 
Thus, the flux relaxation process does not occur in an 
ideal superconductor. 

Let us now check the time dependence of the inter- 
nal field in an ideal superconductor. Putting U c = in 
Ea.pTp. we have 



B(t) = B n 



i=i 



U + t 



(54) 



Eq. ([54| shows that B(t) < 0. Because B(t) is nonneg- 
ative, we must have di = and B(t) = 0. This means 
that the internal field vanishes after the external field is 
reduced to zero. 

Thus, we can conclude that, in an ideal superconduc- 
tor, the flux relaxation process does not occur. The rea- 
son is that the ideal superconductor is free of pinning 
centers, and vortices cannot be pinned down in the bulk 
of the superconductor. 
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B. Inflection point of B(t) — t curve 



inflection point vanishes. 



An earlier study has shown that the B(t)—t curves dis- 
play inflection points ^ The curves have concave shapes 
at short times and then change to convex shapes with 
increasing time. Here I show that this phenomenon does 
not occur when the initial internal field B t reaches a crit- 
ical value. 

For simplicity, let us keep the terms in Eq.(IT5|) to the 
second order B(t) = B e [biw(t) + b 2 w 2 (t)}. (Refill) The 
second derivative of B(t) is 



d 2 B 
dt 2 



= <?(*) 



f 



2b 2 kT 



ln\l + 



U + t 



UpQ 

kT 



(55) 

where g(t) = 2b 2 B e [(r + U + t)/kT]- 2 . 

Letting d 2 B/dt 2 = 0, we obtain the inflection point 
time 



t* 



exp 1 



h Upo 
2b 2 kT kT 



1 



(56) 



Since t* is a time parameter, it must be positive. If 
U > T{exp[l-b 1 /(2b 2 kT) + U p o/(kT)]-l}, then the 



V. CONCLUSION 

In vortex penetration and flux relaxation process, vor- 
tex's activation energy increases with increasing time, 
but its pinning potential does not change. Using this 
idea, one can convert a vortex penetration process with 
a nonzero initial internal field into a process with a zero 
initial internal field by introducing some time parame- 
ters. Similarly, one can also convert a flux relaxation 
process starting with a lower internal field into a pro- 
cess starting with a melting internal field by introducing 
a virtual time interval. This enables us to predict the 
melting internal field from a flux relaxation measurement 
by applying a lower magnetic field to the superconduc- 
tor. Although ideal superconductors do not have pinning 
centers, the vortex penetration into these superconduc- 
tors is still time dependent because of the repulsive force 
of the internal field. But the flux relaxation process does 
not occur in the ideal superconductors. 
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